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Adaptive  Control  of  Uncertain  Hamiltonian  Multi-Input  Multi-Output  Systems: 

With  Application  to  Spacecraft  Control 

Hyungjoo  Yoon  and  Brij  N.  Agrawal 


Abstract — A  novel  adaptive  tracking  control  law  for  nonlinear 
Hamiltonian  multi-input-multi-output  (MIMO )  systems  with  un¬ 
certain  parameters  in  the  actuator  modeling  as  well  as  the  inertia 
and/or  the  Coriolis  and  centrifugal  terms  is  developed.  The  phys¬ 
ical  properties  of  the  Hamiltonian  systems  are  effectively  used  in 
the  control  design  and  the  stability  analysis.  The  number  of  the  pa¬ 
rameter  estimates  is  significantly  lowered  as  compared  to  the  con¬ 
ventional  adaptive  control  methods  which  are  based  on  the  state- 
space  form.  The  developed  control  scheme  is  applied  for  attitude 
control  of  a  spacecraft  with  both  the  inertia  and  the  actuator  un¬ 
certainties,  and  numerical  examples  show  that  the  controller  suc¬ 
cessfully  deals  with  the  unknown  inertia/actuator  parameters. 

Index  Terms — Actuator  uncertainty,  adaptive  control,  Hamil¬ 
tonian  system,  multiple-input-multiple-output  (MIMO)  system 
control,  smooth  projection  algorithm,  spacecraft  attitude  control. 


I.  Introduction 

ADAPTIVE  control  of  multi-input-multi-output  (MIMO) 
dynamic  systems  with  uncertain  parameters  has  been  ex¬ 
tensively  studied  in  the  literature.  Especially,  dynamic  systems 
whose  equations  of  motion  are  described  in  the  general  form 
of  second  order  differential  equations  have  received  significant 
attention  in  the  literature.  Expressing  dynamics  of  systems  in 
the  generic  form,  rather  than  the  state-space  form,  has  several 
advantages.  The  equation  can  be  easily  derived  by  applying  La¬ 
grange’  s  equation,  and  its  form  is  so  general  that  it  can  represent 
various  kinds  of  dynamic  systems,  such  as  a  multilink  robot  ma- 
nipulator[l],  [2]  and  a  spacecraft[3]-[5],  etc.  In  addition,  there 
are  physical  properties,  such  as  energy  conservation,  which  are 
extremely  useful  in  designing  advanced  control  schemes  [1], 
Most  (if  not  all)  of  the  previous  research,  however,  deals  only 
with  uncertainties  in  the  inertia,  centripetal/Corilois,  and  grav¬ 
itational  terms,  assuming  that  an  exact  model  of  the  actuators 
is  available.  This  assumption  is  rarely  satisfied  in  practice  be¬ 
cause  the  actuator  parameters  may  also  have  uncertainties  due 
to  installation  error,  aging  and  wearing  out  of  the  mechanical 
and  electrical  parts,  etc. 

Adaptive  control  with  actuator  uncertainty  does  not  seem  to 
have  received  much  attention  in  the  literature,  even  though  this 
uncertainty  may  result  in  significant  degeneration  of  controller 
performance.  Ge  [6]  has  derived  an  adaptive  control  law  for 
multilink  manipulator  systems  with  uncertainties  in  the  control 
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input  term,  but  the  uncertainty  must  be  in  the  input  scalings,  and 
thus  the  uncertainty  matrix  must  be  diagonal  when  represented 
in  multiplicative  form.  Chang  [7]  has  provided  an  adaptive,  ro¬ 
bust  tracking  control  algorithm  for  nonlinear  MIMO  systems 
which  is  based  on  the  “smooth  projection  algorithm,’’  which  has 
also  been  used  in  [8]  and  [9]  for  adaptive  control  of  SISO  sys¬ 
tems.  More  recently,  using  Chang’s  method  [7],  Yoon  and  Tsio- 
tras  [10]  have  also  provided  an  adaptive  control  scheme  which  is 
applied  to  spacecraft  attitude  tracking  with  uncertain  misalign¬ 
ments/inertia  of  the  actuator  flywheels.  However,  these  previous 
results  [7],  [10]  are  based  on  purely  mathematical  approaches 
and  do  not  exploit  the  useful  physical  properties  of  the  Hamil¬ 
tonian  systems  (In  this  paper,  Hamiltonian  systems  are  referred 
to  as  conservative  systems).  More  significantly,  they  considered 
MIMO  systems  represented  the  state-space  form,  and  thus  they 
need  “over-parameterization”  (more  details  are  in  Section  II). 

In  this  paper,  an  adaptive  tracking  control  algorithm  for 
dynamic  systems  with  parameter  uncertainties  in  the  dynamic 
modeling  as  well  as  the  actuator  modeling  is  developed.  This 
paper  is  organized  as  follows.  Section  II  introduces  the  general 
form  of  second  order  differential  equations  and  the  mathe¬ 
matical  expression  of  the  uncertain  matrices.  In  Section  III, 
adaptive  control  algorithm  is  derived  based  on  the  smooth 
projection  algorithm.  The  proposed  method  is  then  applied  for 
a  spacecraft  attitude  control  problem  with  inertia/actuator  un¬ 
certainties  in  Section  IV.  Finally,  Section  V  provides  numerical 
examples  with  the  spacecraft  to  validate  the  proposed  law,  and 
Section  VI  concludes  this  paper. 

II.  Problem  Formulation 

We  consider  a  MIMO  nonlinear  Hamiltonian  system  repre¬ 
sented  by  the  second-order  differential  equation 

77(q)q+C(q,q)q  +  g(q)  =  Z9(q)u  (1) 

where  q  £  IR"  is  the  generalized  coordinates  vector,  H  £  RraX" 
is  the  (symmetric  positive  definite)  inertia  matrix,  C q  is  a  non¬ 
linear  vector  of  Coriolis  and  centripetal  forces,  and  g  £  R"  is 
the  gravity  vector.  The  generalized  force  is  generated  by  a  con¬ 
trol  input  vector  u  £  R'"  and  the  actuator  matrix  (or  the  control 
influence  matrix)  D  £  R"xm.  For  full  tracking  control  of  q,  it 
is  generally  required  that  n  <  m  and  D  has  full  row  rank. 

Suppose  now  that  the  system  matrices  have  constant  or 
slowly-varying  uncertain  parameters  so  that  they  can  be  as¬ 
sumed  to  be  constant.  Let  us  express  these  matrices  as 

H( q,  0S)  =  Hn{ q)  +  Ha( q,  0S)  (2a) 

C(q,  q,  0S)  =  Cn(  q,  q)  +  CA(q,  q,  0S)  (2b) 

g(q,  0S)  =  g,l(q)  +  gA(q,  0S)  (2c) 

D( q,  0a)  =  Dn( q)  +  DA(q,  0O)  (2d) 
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where  the  matrices  with  a  superscript  of  “n”  are  with  their 
known  nominal  values,  Oa  £  Rp  is  a  vector  of  unknown 
bounded  constant  uncertainties  in  the  actuator  matrix  D , 
and  0,  £  R'1  is  in  the  other  system  matrices/vector,  H, 
C,  and  g.  We  also  assume  that,  by  proper  definition  of  the 
unknown  parameters  Qs.  and  @a,  the  uncertain  matrices 
iTA(q,0,),CA(q,q,0,),gA(q,0s)  and  DA(q,0a)  are 
linearly  dependent  on  0S,  and  0a,  respectively. 

The  problem  we  study  in  this  paper  is  formulated  as  follows. 

1 )  Adaptive  Tracking  Control  of  Uncertain  MIMO  System: 
Consider  the  MIMO  Hamiltonian  System  ( 1 )  with  system  ma¬ 
trices  with  uncertainties  as  expressed  in  (2).  For  given  reference 
trajectory  q d(t),  qd(t)  and  qd(t)  which  are  bounded,  design  a 
control  algorithm  which  makes  the  tracking  error  q  =  q  —  q^ 
globally  asymptotically  stable  at  q  =  0. 

Slotine’s  adaptive  controlled  1]  can  handle  only  the  special 
case  with  DA  =  0.  On  the  other  hand,  Chang’s  method  [7] 
needs  to  convert  the  equation  of  motion  (1)  into  a  state-space 
form,  that  is 

q  0nxri  I n  q 

q  0 nXn  0 nXn _  q_ 

+  °"x"  (-H~lCii-H-1e  +  H-1Da)  (3) 

J-n 

in  which  the  uncertainties  in  the  inertia  matrix  H  are  multiplied 
with  the  uncertainties  in  the  other  terms.  Since  Chang’s  method 
also  uses  the  linear  dependency  of  uncertainties,  it  needs  to  es¬ 
timate  the  combination  of  uncertain  parameters  in  H  and  those 
in  the  others  (“over-parameterization”).  Therefore,  the  number 
of  parameter  estimates  would  significantly  increase  when  H  has 
multiple  uncertain  parameters.  In  the  following  section,  we  de¬ 
rive  a  novel  algorithm  which  remedies  these  drawbacks,  that  is 
an  adaptive  tracking  controller  for  a  case  of  DA  f  0  without 
over-parameterization. 

III.  Adaptive  Control  Law 

The  first  part  of  the  derivation  of  the  adaptive  control  law 
follows  the  standard  design  procedure  for  Hamiltonian  systems 
in  [1]  and  [4].  Let  us  define  a  measure  of  tracking  error  s  = 
q  +  Aq  =  q  —  q?,  and  the  reference  velocity  q,,  =  qd  —  Aq, 
where  the  matrix  —A  is  assumed  to  be  Hurwitz.  Let  0*  be  the 
parameter  estimate  vector  and  let  0*  =  0*  —  0*  be  a  parameter 
estimate  error  vector,  when  *  is  a  or  s. 

As  suggested  by  Slotine  et  al.  [1]  from  a  physical  insight  that 
q7//q  is  the  system’s  kinetic  energy,  the  following  Lyapunov 
function  candidate  is  defined  as: 

V(t)  =  \  [st7Ts  +  ©Jr;1©,  +  ©Jr;1©,]  (4) 

where  Ta,  Fs  are  positive  definite  weighting  matrices.  Differen¬ 
tiating  V ( t )  with  respect  to  time  and  using  the  skew-symmetry 
property  of  the  matrix  H  —  2 C  to  replace  the  term  (1  /2)s77/s 
with  sTCs,  one  can  have  the  following  expression: 

V—  stHs  +  ± stHs  +  ©Jr;1^),  +  ©Jr;1^, 

=  sT  (Hs  +  Cs)  +  ©Jr;1^,  +  ©Jr;1^ 
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=  sT(JDu-jTqr-Cq?.-g)+0rr-14+efrj14)s.(5) 

When  the  actuator  matrix  D  is  exactly  known  (i.e.,  DA  = 
0)  and  has  full  row  rank  (i.e.,  rank(D)  =  n),  one  can  easily 
design  an  adaptive  control  law  using  the  methods  proposed  in 
the  previous  work  [1].  However,  since  D  is  assumed  to  contain 
unknown  parameters  as  well  as  H,  C  and  g,  a  novel  control  law 
is  proposed  as  follows: 

(Dn  +  Da)  u  =  [{Hn  +  Ha)  qr 

+  (Cn  +  C'A)qr  +  (g"  +  gA)  —  Kds'j  (6) 

where  the  matrices/vector  with  a  ‘hat’  symbol  are  constructed 
using  the  parameter  estimates  0a  and  0S  instead  of  the  (un¬ 
known)  actual  parameters.  Kd  is  a  gain  matrix  which  is  a  posi¬ 
tive  definite.  When  the  matrix  Dn  +  D A  is  assumed  to  have  full 
row  rank,  the  (weighted)  minimum  norm  solution  is  given  by 

u  =  ( Dn  +  Z)A)t  [{Hn  +  Ha)  qr 

+  (Cn  +  CA) qr  +  (g"  +  gA)  -  Kd s)  (7) 

where  (-)t  denotes  the  (weighted)  pseudo-inverse  of  a  matrix 

[5],  [10]. 

The  control  law  (7)  leads  to 

V  =  sr[{(Dr!  +  Da)  -  DA)}u  -  (Hn  +  Ha) q?. 

-  (C"+cA)qr  -  (g”+gA)]+0jr-14)a  +  ©^rj1^ 
=  sT[(iT"  +  Ha)  qr  +  (Cn  +  CA)  qr  +  (gn  +  gA) 

-  Kds  -  Dau  -  (Hn  +  Ha) qr  -  (Cn  +  CA) qr 

-  (gil  +  gA)]  +  ejr -14  +  ©^rj1^, 

=  sT  Ha qr  +  CAqr  +  gA  -  Kd s  -  DA u 

+  QTaT-^a  + 

=  -  sTKds  +  sT(HAqr  +  CAqr  +  gA) 

-  sT  dau  +  ©Jr -10a  +  ©Jr;1©,.  (8) 

Since  the  uncertainty  matrices/vector  are  assumed  to  depend 
linearly  on  0’s,  one  can  define  known  regressor  matrices  (in 
fact,  row  vectors)  Ys  =  Y,(q,  q, q,., q,., s)  £  R1X9  and  Ya  = 
Y„(q,  q,  qr,  qr.  s,  u)  £  Rlxp,  such  that 

sT(HAqr  +  CAqr  +  gA)  =  Y,0,  (9) 

and 

-stDau  =  Y„0„.  (10) 

Notice  that  the  definition  of  the  regressor  matrix  Ys  in  (9)  is 
slightly  different  from  that  in  the  previous  works  by  Slotine  [1], 
[3],  in  which  the  regressor  matrix  was  defined  as  ( HAqr  + 
CAq,,  +  gA)  =  Ys<ds.  By  including  s  in  the  regressor’s  defi¬ 
nition,  the  size  of  the  regressor  Ys  becomes  smaller  than  that  of 
Slotine’s  works  where  the  regressor  has  a  size  of  n  X  q. 

The  time  derivative  of  V  then  becomes 

v  =  -sT7rds+Ys0s+Ya0a+0jr;14)s+0jr-14)a  an 
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and  taking  the  adaptation  laws  of  the  parameter  estimates  to  be 

4s  =  -rsysT  (12) 

and 

4,  =  -rai;T  d3) 

then  yields  V(t)  =  —sTKds  <  0.  Using  standard  arguments  in 
[  1],  [10]  which  use  Barbalat’s  lemma,  one  can  easily  show  that 
V  — >  0  and  thus  s  — >  0  and  q  — >  0  as  t  — *  oo. 

A.  Smooth  Projection  Algorithm 

We  previously  used  an  assumption  that  the  matrix  Dn  +  DA 
has  full  row  rank  in  deriving  the  adaptive  control  laws  (7),  (12), 
and  (13).  In  general,  for  the  full  tracking  control,  the  nominal 
matrix  Dn  in  general  has  full  row  rank.  However,  a  drift  of  the 
parameter  estimates  0a,  governed  by  an  update  law  (13),  can 
result  in  Dn  +  T)A  losing  rank.  We  will  refer  to  this  situation 
as  a  “singularity”  of  the  steering  law  due  to  the  adaptation.  This 
singularity  hinders  the  use  of  the  derived  control  laws,  so  they 
need  to  be  modified. 

If  the  nominal  matrix  Dn  has  full  row  rank,  and  the  true  value 
of  the  parameter  uncertainty  Qa  is  bounded  by  a  sufficiently 
small  number,  and  the  parameter  estimate  0a  is  also  kept  small, 
then  the  matrix  Dn  +  DA  will  also  have  full  row  rank.  To  this 
end,  we  define  the  following  two  convex  sets: 

^  {©a  e  R^l  ||©o||2  </?}  (14) 

^©o  =  {BaGRp|||0a||2</3  +  5}  (15) 

where  /3  >  0  and  S  >  0  are  known  constants.  Notice  that  0©a  C 
ClQa .  We  make  the  following  three  assumptions. 

•  Assumption  1.  The  nominal  value  D"  has  full  row  rank  of 
n. 

•  Assumption  2.  The  actual  value  @a  belongs  to  the  set  . 

•  Assumption  3.  If  0a  G  ,  then  D"  +  DA  is  non¬ 
singular. 

These  assumptions  allow  us  to  modify  the  adaptation  law  (13) 
by  using  the  “smooth  projection  algorithm”  as  follows:1 

©a  =  Proj(0a,$a)  (16) 

where 

=  - raY( J  (17) 

'The  adaptation  law  is,  in  fact,  only  Lipschitz  continuous,  not  continuously 
differentiable.  The  use  of  the  term  "smooth”  is  a  slight  misnomer  in  this  context, 
but  we  use  it  here  in  accordance  to  prior  usage  in  the  literature.  It  should  be 
noted  that  a  new  parameter  projection  operator  which  is  Cn  has  been  recently 
introduced  in  [11], 


and  shown  in  (18)  at  the  bottom  of  the  page.  This  adapta¬ 
tion  law  is  identical  to  (13)  in  cases  (i)  and  (ii),  and  switches 
smoothly  to  a  new  expression  in  case  (iii).  The  projection  op¬ 
erator  Proj(0a,  <!>„)  is  locally  Lipschitz  in  (0a,  <!>„).  thus  the 
system  has  a  unique  solution  defined  for  some  time  interval 
[0 ,T),  T  >  0. 

Proposition  1:  Under  Assumptions  1-3,  the  control  law  (7) 
along  with  the  adaptation  laws  (12)  and  (16)  yields 

V  <  -s TKds  <  0  (19) 

and 

©a(f  =  0)  G  flea  =>■  Qa(t)  G  H0o,  Vi>0.  (20) 

Proof:  The  proof  is  straightforward  and  therefore  omitted 
here.  It  is  similar  with  the  proof  in  the  author’s  previous  work 
[10].  ■ 

From  Proposition  1,  one  can  conclude  that,  using  the  feed¬ 
back  control  law  (7)  and  the  adaptation  laws  (12)  and  (16), 
q  — >  0  as  t  — ■»  oo  and  ( Dn  +  DA )  will  not  lose  rank,  if  we 
choose  the  initial  parameter  guess  0a(O)  inside  the  set  floa. 
For  instance,  we  may  take  0a(O)  =  0. 

It  is  also  worth  mentioning  that  the  proposed  adaptation  law 
(16)  has  the  additional  benefit  of  keeping  the  parameter  esti¬ 
mates  from  “bursting,”  which  may  happen  when  the  persistency 
of  excitation  condition  does  not  hold  [12], 

IV.  Application  to  Spacecraft  Attitude  Control 
A.  Equations  of  Motion 

In  this  section,  applying  the  proposed  adaptive  control 
scheme,  we  design  an  adaptive  attitude  tracking  control  law 
for  a  spacecraft  which  has  uncertainties  both  in  its  inertia  and 
actuator  modeling.  (Since  the  system  remains  Hamiltonian 
in  the  presence  of  uncertainties,  non-adaptive  passivity-based 
methods  [13],  [14]  may  be  applied  to  control  the  spacecraft 
attitude,  but  they  are  regulating  control  laws  which  can  be  used 
only  when  q d(t)  is  constant.) 

A  cluster  of  variable-speed  control  moment  gyros  (VSCMGs) 
with  N  flywheels  is  used  for  the  torque  actuator.  While  a  con¬ 
ventional  control  moment  gyro  (CMG)  keeps  its  flywheel  spin¬ 
ning  at  a  constant  rate,  a  VSCMG — as  its  name  implies — is 
essentially  a  single-gimbal  CMG  with  the  flywheel  allowed  to 
have  variable  speed.  (See  [15]  and  [16]  for  more  details  and  ap¬ 
plications  of  VSCMGs.) 

Fig.  1(a)  shows  a  spacecraft  with  the  ith  VSCMG,  where  g, 
is  (body-fixed)  gimbal  axis,  s,  is  spin  axis,  and  t ;  =  g;  X  s7;  is 
transverse  torque  axis.  The  equations  of  motion  of  a  spacecraft 
with  VSCMGs  are  complicated  as  shown  in  the  aforementioned 


Proj(0a,<t>a) 


(l|eall2-/3)Ce, 

«ll©<.ll2 


if  (i)  ||0a||2  </3,  or  (ii)  ||0a||2  >  8  and  ^  Qa  <  0, 
if  (iii)  ||0a||2  >  8  and  $a  ©a  >  0. 


(18) 
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(a)  Spacecraft  body  and  the  z'-th  VSCMG  frames  (b)  The  Pyramid  Configuration 


Fig.  1.  Spacecraft  and  VSCMGs. 


references,  but  under  assumptions  which  are  standard  in  the  lit¬ 
erature  [10],  they  can  be  simplified  as  follows: 

Jut  -  [hx]o;  +  Qu  =  0  (21) 

and 

h  =  Jlj  +  AsIwCl  (22) 

where  Q  =  [Qcm.g,Qrw]  £  R3x2JV  and  where  Qcmg{ 7,0)  = 
Atlwfld  £  R3xiv,  Qrwi'Y)  =  ASIW  £  [R3xiv,  and  the  control 
input  of  this  system  is  u  =  [u\, . . . ,  U2n]T  =  [yLft  ]T  £ 
R2A  .  J  is  the  total  moment  of  inertia  of  the  spacecraft  which 
is  assumed  to  be  constant,2  cj  is  the  body  rate  vector  of  the 
spacecraft,  h  is  the  total  angular  momentum  of  the  spacecraft, 
7  =  [71, . . .  ,7jv]T  £  RA  and  Cl  =  [S2i, . . . ,  Qn]t  £  RA  are 
vectors  of  gimbal  angles  and  flywheel  spinning  speeds,  respec¬ 
tively,  and  Iw  is  a  diagonal  matrix  with  the  inertias  of  VSCMGs 
flywheels.  The  skew-symmetric  matrix  [vx],  for  v  £  R3,  rep¬ 
resents  the  cross  product  operation.  The  matrices  A*  £  R3xAr 
have  as  columns  the  gimbal  (g.,),  spin  (sA,  and  transverse  (t.,) 
directional  unit  vectors  expressed  in  the  body-frame,  where  *  is 
g,  s  or  t.  These  matrices  depend  on  the  gimbal  angles  as  follows: 


and  Ir  is  the  r  X  r  identity  matrix. 

As  suggested  in  [1]  and  [5],  we  combine  the  kinetic  (21)  and 
the  kinematic  (26)  into  one  second-order  system  as  follows: 

H(i 7)0  +  C(er,  0)0  =  D{a)  u  (28) 


where 

H(0)=G-tJG~ 1  (29) 

C(0,&)  =  -  G_TJG_1GG_1  -  G-T[(i?(a)h7)x]G-1  (30) 

D(a)=-G~TQ.  (31) 

R{0)  is  a  rotational  matrix  from  the  inertial  frame  to  the  body 
frame,  and  h/  is  the  total  angular  momentum  of  a  spacecraft 
expressed  in  the  inertial  frame  which  is  conserved  to  be  constant 
if  there  is  no  external  torque  applied  to  the  spacecraft.  Therefore, 
h  =  R(a)  h7. 

Notice  that  the  equation  of  motion  (28)  has  the  form  of  (1) 
with  the  gravitational  term  g  =  0.  Moreover,  it  can  be  easily 
shown  that  the  matrix  H  —  2 G  is  skew-symmetric[  1], 

B.  Adaptive  Attitude  Tracking  Control 

Suppose  that  there  are  uncertainties  in  D  as  well  as  J.  We 
assume  that  the  exact  values  of  the  initial  axis  directions  of 
VSCMGs  actuator  at  7  =  0  are  unknown.  This  can  happen  when 
the  VSCMGs  are  installed  with  small  misalignments  and/or  the 
measure  of  gimbal  angles  has  constant  unknown  bias.  In  addi¬ 
tion,  h7  is  also  unknown  constant  not  only  because  of  uncertain 
J  but  also  because  of  uncertain  As. 

For  most  cases  the  effect  of  axes  uncertainties  on  the  overall 
system  performance  may  be  not  significant.  However,  when  fly¬ 
wheels  spin  at  high  speeds  as  in  a  case  of  integrated  power  and 
attitude  control  system  (IPACS)  [5],  [18],  [19],  even  small  mis¬ 
alignments  of  the  flywheel  axes  can  be  detrimental. 

The  uncertain  parameters  in  H  and  C  can  be  defined  as  fol¬ 
lows: 


Ag  —  Ago  (23) 

As  =  As0  [cos  7]^  +  Ato  [sin  if  (24) 

At  =  Ai0[cos7]d  -  As0[sin7]d  (25) 

where  the  A*o’s  denote  the  values  of  A*  at  7  =  0.  The  symbol 
x7  denotes  the  diagonal  matrix  with  elements  the  components 
of  the  vector  x,  and  cos7  =  [cos  71, . . . ,  cos7jv]t  and  sin7  = 
[sin7i, . . . ,  sin77v]r. 

The  modified  Rodrigues  parameters  (MRPs)  [17]  are  chosen 
to  describe  the  attitude  kinematics  of  the  spacecraft.3  The  kine¬ 
matics  in  terms  of  the  MRPs  is  given  by 

0  =  G{o)  u  (26) 


where 


G{tr)  =  *  (l3  +  [0X]  +  00T  -  *(1  +  0 


(27) 


2In  fact,  J  is  a  function  of  7  but  the  dependence  is  weak  in  general. 

3We  hasten  to  point  out  that  the  use  of  the  MRPs  to  describe  the  kinematics  is 
done  without  loss  of  generality.  Any  other  suitable  kinematic  description  could 
have  been  used  with  the  conclusions  of  the  paper  remaining  essentially  the  same. 


©s  —  [Ajn,  Aj22,  Aj33,  Aji2,  Aji3,  Aj23,  A/li,  A/l2,  a hs]T 
£  R9  (32) 

where  Aj's  are  the  elements  of  JA  =  J  —  Jn  and  J"  is  the 
nominal  value  of  the  actual  inertia  matrix  J .  Similarly,  A/i’s 
are  elements  of  hJA  =  h7  —  h7ra .  The  uncertain  parameters  in 
D  are  defined  as 


1 )6N 


(33) 


where 


©t,i  =  ti,0  —  t”0, 


(-) , :  —  s,:n  s 


i,0  sifli 


i  =  l,...,N.  (34) 


and  t.,;.o  and  s;.o  are  actual  value  of  t,  and  s7;  at  7^  =  0,  re¬ 
spectively,  and  t"0  and  s”0  are  their  nominal  values.  The  total 
number  of  parameter  estimates  is  then  6 N  +  9.  If  one  uses  the 
methods  in  the  previous  works  [7],  [10],  then  the  number  of  the 
estimates  will  be  as  much  as  6  X  (6 N  +  3). 

Exact  mathematical  expressions  of  the  regressor  vectors  Ys 
defined  in  (9)  can  be  easily  obtained  using  symbolic  math  pack¬ 
ages,  and  can  be  constructed  from  the  measurements  of  0,  0, 
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and  the  desired  trajectories  a,/-  a,i,  ij,i.  The  regressor  vector  Ya 
defined  in  (10)  can  be  obtained  in  the  same  way,  but  it  is  also 
possible  to  derive  its  mathematical  expression  by  manipulating 
the  matrices  as  follows: 

Ya  =  stG-t[(Mu)1I.3,  . . . ,  (Mu)2JVI3]  (35) 


which  is  obtained  by  adding/subtracting  20%  of  the  nominal 
values.  The  (unknown)  actual  value  of  Qs  defined  in  (32)  is  then 

0S  ~  [3000,  -1300,  2400,  -600,  -200, 

-400,  -89.77,  141.43,  -162.77]T. 


where 


[cos7]d 

.s 

.‘53. 

"  Lwnd  0 " 

—  [sin  7] d 

[cos7]d 

0  Iw 

(36) 


and  (x)i  is  the  7th  element  of  a  vector  x.  Then  using  the  devel¬ 
oped  control  scheme  in  Section  III,  one  can  design  an  adaptive 
attitude  tracking  control  law  for  the  spacecraft. 


V.  Numerical  Examples 

Numerical  examples  for  a  satellite  with  a  VSCMGs  cluster 
are  provided  in  this  section  to  test  the  proposed  adaptive  con¬ 
trol  algorithm.  In  addition,  the  null  motion  technique  presented 
in  [19]  is  also  applied  along  with  the  proposed  control  law  in 
order  to  avoid  the  geometric  singularity  of  VSCMGs.  Since  the 
control  input  u  in  (7)  is  replaced  by  u+unun,  where  the  null  mo¬ 
tion  unun  satisfies  ( Dn  +  Z)A)unun  =  0,  the  stability  analysis 
of  the  adaptive  control  law  still  holds.  A  standard  four-VSCMG 
pyramid  configuration  (N  =  4)  is  utilized  as  shown  in  Fig.  1(b). 

When  the  skew  angle  8  in  Fig.  1(b)  is  chosen  so  that 
cos  8  =  l/\/3,  the  nominal  values  of  the  axis  directions  at 
7  =  [0, 0, 0,  Of  are 


0  -1 

0  r 

1  0 

-1  0 

(37) 

0  0 

0  0_ 

-0.5774 

0 

0.5774 

0 

0 

-0.5774 

0 

0.5774 

.  (38) 

0.8165 

0.8165 

0.8165 

0.8165 

The  (unknown)  actual  axis  directions  at  7  =  0  used  in  the 
present  example  are  assumed  as 


-0.0072 

-0.9999 

-0.0072 

0.9999 

Aso  = 

0.9999 

0.0143 

-0.9999 

0.0057 

-0.0071 

0.0071 

-0.0142 

-0.0117 

' -0.5657 

-0.0024 

0.5831 

0.0061" 

Ato= 

0.0018 

-0.5715 

-0.0157 

0.5868 

0.8246 

0.8206 

0.8123 

0.8097 

which  are  obtained  by  rotating  each  nominal  gimbal  frame 
(g;,  s,,  t,j  with  1  degree  about  arbitrary  direction.  With  these 
values,  ||0a||2  —  0.0015.  The  nominal  value  of  the  spacecraft 
inertia  matrix  is 


"  15000 

3000 

-1000" 

Jn  = 

3000 

6500 

2000 

-1000 

2000 

12000 

(unknown)  actual  inertia  matrix  is 

18000 

2400 

-1200" 

J  = 

2400 

5200 

1600 

-1200 

1600 

14400 

kg  m2 


kg  m2 


(41) 


(42) 


The  remained  parameters  used  for  the  simulations  are 
set  as  follows:  the  initial  attitude  is  aligned  with  the  in¬ 
ertia  frame  at  rest  (<r(0)  =  o)(0)  =  ot(0)  =  [0, 0, 0]T), 
the  initial  gimbal  angle  is  also  zero  at  rest  (7  =  7  = 
[0, 0, 0,0]T),  the  moment  of  inertia  of  the  flywheels  are 
Iw  =  diag{2.0, 2.0, 2.0, 2.0}  kgm2,  the  initial  wheel  spin 
speed  is  0(0)  =  104  X  [2.5, 3.5, 3.5, 3. 0]r  rpm,  and  the 
controller  parameters  are  chosen  as  Kd  =  1 0'T',,  A  =  I3, 
/ 3  =  0.01,  and  6  =  0.01.  (Notice  that  Assumption  2  is  satisfied 
with  the  given  value  of  /3 .) 

Notice  that  the  initial  wheel  speeds  of  the  VSCMGs  are  set  to 
25  000  ~  35  000  r/min,  which  are  an  order  of  magnitude  larger 
than  the  speed  of  conventional  CMGs,  since  the  flywheels  of 
VSCMGs  used  for  IPACS  in  general  need  to  spin  very  fast  so 
that  they  are  competitive  against  traditional  chemical  batteries. 
According  to  [20],  even  a  higher  speed  than  these  vales  is  im- 
plementable,  at  least  in  a  laboratory. 

The  reference  trajectory  is  chosen  so  that  the  ini¬ 
tial  reference  attitude  is  aligned  with  the  body  frame 
which  is  also  aligned  with  the  inertial  frame,  and  the 
angular  velocity  of  the  reference  attitude  is  chosen  as 
ujd(t)=  0.04(sin  27rf/400,  sin  27rf/300,  sin  27rf/200)rrad/s. 

First,  in  order  to  show  the  effect  of  the  misalignment  of  the 
axes  of  the  VSCMG  cluster  and  the  uncertain  inertia  matrix,  a 
simulation  without  adaptation  was  performed.  (This  nonadap- 
tive  control  law  can  be  regarded  as  a  combination  of  a  feed¬ 
forward  term  and  a  simple  proportional-derivative  (PD)  control 
term  [1].)  Fig.  2(a)  shows  the  attitude  tracking  error  (expressed 
with  “3-2-1”  Euler  angles)4  under  control  law  with  the  adap¬ 
tation  gains  Ta  and  Fs.  set  to  zero  matrices.  Since  the  flywheel 
speeds  are  very  fast,  there  is  large  attitude  tracking  error  without 
adaptation. 

Next,  another  simulation  was  run  with  adaptation  of  the  actu¬ 
ator  uncertainty  0a  only.  The  adaptation  gains  are  set  to  F„  = 
10  IqN  and  Fs  =  0,  and  the  resulting  attitude  error  is  shown  in 
Fig.  2(b).  There  is  significant  performance  improvement  with 
the  adaptation  of  0a  only,  but  tracking  error  with  a  magnitude 
of  about  0.2  degree  remains.  On  the  other  hand.  Fig.  2(c)  shows 
the  attitude  tracking  error  with  adaptation  of  0S  only.  The  adap¬ 
tation  gains  are  Ta  =  0  and  Ts  =  diag(10'  Ig,  105  I3).  In  fact, 
the  control  law  in  this  scenario  is  almost  identical  with  Slotine’s 
method  [1],  The  attitude  error  is  significantly  attenuated  using 
the  adaptive  controller,  but  there  are  again  residual  tracking  er¬ 
rors  with  a  magnitude  of  about  0. 1  degree. 

Finally,  a  simulation  is  performed  with  adaptation  of  both  I), 
and  Ts,  that  is  Ta  =  10  Igjv  and  L,  =  diag(107  Iq,  105  I3). 
Fig.  2(d)  shows  the  tracking  performance  is  significantly  im¬ 
proved  upon  Slotine’s  control  law.  Fig.  3(a)  shows  the  time 
history  of  ||0a||2.  It  is  confirmed  that  ||0a||2  does  not  drift 

4The  use  of  Euler  angles  in  the  figures  in  done  solely  for  the  convenience  of 
the  reader  who  may  not  be  familiar  with  the  MRPs. 
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Attitude  Error  Attitude  Error 


Attitude  Error  Attitude  Error 


Fig.  2.  Attitude  tracking  errors. 


(a)  Estimates  of  Qa  (square  of  norm) 


(b)  Estimates  of  A J 


(c)  Estimates  of  Ah 


Fig.  3.  Parameter  estimates. 


more  than  (3  +  S  =  0.02  owing  to  the  smooth  projection  al¬ 
gorithm.  As  a  result,  the  steering  law  (7)  remains  well-defined. 
Fig.  3(b)  and  (c)  show  the  time  history  of  the  parameter  esti¬ 
mate  Qs.  In  Fig.  3(c)  the  bold  horizontal  dotted  lines  denote  the 
actual  values  of  the  components  of  lr,A.  Some  of  the  estimates 
approach  their  actual  values,  but  they  do  not  exactly  converge 
to  these  values  due  to  the  lack  of  persistency  of  excitation.  Note 
that,  in  general — and  for  a  linear  system — convergent  estima¬ 
tion  of  m  parameters  requires  at  least  m/2  sinusoids  in  the  ref¬ 
erence  signal.  For  the  nonlinear  case  such  simple  relation  may 
not  be  valid  [1],  In  this  example,  the  number  of  the  parame¬ 
ters  to  be  estimated  is  6N  +  9,  while  the  reference  signal  has 
only  three  sinusoids,  thus  the  persistency  condition  is  not  sat¬ 
isfied.  However,  the  tracking  error  will  eventually  converge  to 
zero  nevertheless. 

VI.  Conclusion 

In  this  paper,  we  proposed  an  adaptive  tracking  control  law 
for  a  nonlinear  Hamiltonian  MIMO  dynamic  system.  The  pro¬ 
posed  control  scheme  has  several  significant  improvements  over 


the  previous  works  in  the  literature.  First,  the  proposed  method 
can  deal  with  uncertainties  in  the  actuator  terms.  Second,  the 
proposed  method  exploits  the  physical  properties  of  the  Hamil¬ 
tonian  systems  and  so  the  designed  law  is  more  compact  than 
those  in  previous  mathematics-based  methods.  Finally,  it  does 
not  need  over-parameterization  to  deal  with  uncertainties  both 
in  the  inertia  and  the  actuator  modeling  at  the  same  time,  while 
the  previous  methods  do. 

The  developed  adaptive  algorithm  is  shown  to  signifi¬ 
cantly  improve  the  tracking  performance  in  the  application 
to  the  spacecraft  attitude  control,  but  it  still  has  room  for 
improvement.  For  instance,  while  only  three  parameters  are 
generally  needed  to  express  a  misalignment  of  axis  frame  for 
one  VSCMG,  a  total  of  6  parameters  are  used  in  this  method. 
Development  of  methods  to  reduce  the  number  of  estimated 
parameters  would  be  extremely  beneficial.  In  addition,  robust¬ 
ness  of  the  proposed  control  law  against  unknown  external 
disturbances  is  also  suggested  for  future  study. 
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